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1. THE MA-088 HGD DATA 


Of the two 5 84 A channels in the MA-088, namely #1 and #3, channel 
#1 appeared to provide data with erratic count rates and undue susceptibility 
to day glow and solar contamination possibly because of filter fatigue or 
failure. Channel #3 data appear normal and of high quality. For this 
reason only data from this last channel was analyzed and used for detailed 
comparison with theory. All MA-088 data have been binned into 30-s 
interyals , corresponding to the time necessary for the LOS to traverse 
at most one field width. During each such interval we separately treated 
data taken when, the gas absorption cell was full and empty, and ignored 
data taken while the cell pressure was changing. In order to eliminate 
spikes of extraneous noise and telemetry dropouts from the data, we 
subjected each bin of 0.1-s data samples to a filtering process in which 
the mean and standard deviation of each set of samples was calculated, and 
then all samples more than three standard deviations away from the mean 
were discarded. •• The process was repeated until no more, samples were 
discarded. The amount of data thereby eliminated was less than five per 
cent and was consistent with Poisson statistics. The resulting count rates 
for the various bins ranged between 2 and 60 counts s ^ at night and 12 
to 1880 counts s ^ in the daytime . The standard deviations obtained from 
the -statistics of the individual 0.1-s data samples were consistent with 
those expected from Poisson statistics. 

We corrected the observed count rates for atmospheric absorption of 

O 

584 A flux along the lines of sight by using the atmospheric models of 
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Jacchia (19.71) and .the 584 A absorption cross-sections of Hall, Schweizer, 
and Hinteregger (1965). for N 2 , 0 and 0^. Our calculations were based on 
the spacecraft position and LOS orientation at the mid-point of each bin. 
The solar activity parameters required as inputs to the Jacchia models 
were taken from Lincoln (1975) , The effect of these corrections, was to 
increase the observed flux levels by 10 to 30 per cent. 

O 

It is also necessary to allow for resonant absorption of 584 A flux 
14 -2 

by the 1-2 x 10 atoms cm (Jacchia 1971) of atmospheric helium remaining 
above the spacecraft. This correction is model-dependent, for it depends 

O 

on the spectral profile of the 584 A line scattered from the ISM helium. 
Therefore, we chose to incorporate this correction into the models used- 
for predicting the scattered line. 


2. APPLICATIONS Of THE HGD DATA 

2."1 .The Nighttime 584 A Emission. 

To make certain that our conclusions about flux scatterd from the 

' O 

local interstellar medium at night were not affected by daytime 584 A 
radiation, we restricted the analysis in this case to data taken when the 
Sun was at least 30° below the horizon. Since the Earth’s limb is depressed 
15° below the horizontal when viewed from 220 km altitude, this constraint 
corresponded to a solar zenith angle of at least 135°. To minimize the 

O 

affect on our analysis of errors in calculating absorption of 584 A flux 
■by the atmosphere above 220 km, we chose to analyze only those data for 
which the zenith angle of the LOS was less than 60°. Finally, we excluded 
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from analysis data taken whenuthe spacecraft was in the South Atlantic 
Anomaly-. After these three constraints were imposed, 53 minutes of data 
• remained for study . 

We investigated the possibility of other terrestrial phenomena 
.affecting the data by correlating the binned, absorption-corrected data 
with LOS zenith angle, solar zenith angle, LOS-Sun angle, local strength 
of the geomagnetic field, magnetic ..latitude, and component of spacecraft 
velocity along the LOS . The only correlation we found at night was with 
magnetic latitude . There is a pronounced increase in count rates south 
of magnetic latitude -50°, and a noticeable peak within 10° of the magnetic 
equator. These increased count rates may be due to localized airglow such 
as that reported by Carruthers and Page (1972), or to some other cause. 

We decided to restrict the data under analysis further, by excluding data 
taken south of magnetic latitude -50° or within 10° of the magnetic 
equator. After this additional exclusion there remained 42 minutes of 
data. 

We encountered one more problem in reducing the data. -As discussed 

O 

by Schmidtke (1976), the total solar 584 A flux varies from, day to day 
by as much as thirty per cent. This short-term variability appears to be 
sufficiently unsystematic that we are reluctant to interpolate between 
Schmidtke' s published data points, which are spaced at least 24 hours 

O 

apart. We are not aware of any studies of solar 584 A variability on 
shorter time scales. 

The 42 minutes of Apollo-Soyuz data we have selected for analysis 
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are composed of eight separate intervals, 1-10 minutes long, spaced 1-24 

O 

hours apart. To allow for possible solar 584 A flux variations from 
interval to interval,, we decided' to allow the assumed value of that flux 
to vary, from interval to interval, while searching for model fits. 

The by now conventionally accepted modified cold model (MCM) of ISM 
flow through the solar system was used to attempt to interpret these data 
Via model-fitting. The cold model phase-space density function derives 
from the assumption that the distant ISM has a temperature of zero; that 
is, that all atoms move in the same direction with the same speed. Given 
a point near the Sun, not located on the axis of symmetry, it can be 
shown that only two orbits with the zero— tenperature initial conditions 
cross that point. 

If the selected point is directly upstream of the Sun only one orbit 
crosses it, but in the zero-temperature approximation every orbit crosses 
the downstream symmetry axis. This latter condition leads to an unphysical 
density spike at the downstream axis, and the cold model thereby falsely 
predicts a bright, localized EUV flux scattered from that region. 

Feldman, Goldstein and Sherb (1972) and Blum, Pfleiderer and Wulf-Mathies 
(1975) have published formulae for altering the number density function 
close to the downstream axis to take first-order account of non-zero 
atomic thermal speeds. With the addition of such formulae, the cold model 
becomes the modified cold model. 

In this model, the three-dimensional integral over gas velocity 
collapses to a two-term sum, and Dirac delta functions allow the integrations 
on frequency to be performed immediately. The evaluation of the flux I in 
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this model is thereby reduced to a one-dimensional numerical integration 
along the line of sight’, which can be accomplished quickly on even a 
small computer. 

We have used this modified cold model to attempt to explain the HGD 
observations.' No reasonable fits to the data were found. This is 
probably due to the fact that the cold model is valid only when the thermal 
gas speeds are small compared to the bulk speed. The ASTP HGD instrument 
is highly sensitive, however, to the dynamical state of- the gas. Thus 
the inability of the cold model to explain the data must be taken to 
indicate that the instrument is far too sophisticated for such a simple 
model and to make any progress a better model must be used. 

We have developed this model and have applied it to the HGD observations. 
A brief description of this model (which we call the hot model) is given 
in paragraph 3 of this report and a more detailed analysis of it may be 
found in paper 6.3 in-.the appendix. 

The hot model predictions were again compared to the HGD data in a 
manner similar to that used unsuccessfully with the cold .model. Details 
of this analysis may be found in paper 6.1 in the appendix to this report. 

In brief, we evaluated the hot model predictions for the physical 
parameters corresponding to the mid-points of the various HGD data bins. 

In all cases, the cell-full count rates and predictions *were substantially 
less than the corresponding cell-empty values . Almost all cell-empty flux 
predictions . increased with temperatures as did most cell-full predictions. 
During one interval of data collection, the cell-full predictions were 
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extremely low, between 0.001 and 0.06 ray leigh, ' while the actual count 
rates detected by the instrument were 2 or 3 per second. If these count 

O 

rates were caused by helium 584 A flux, they correspond to fluxes of 
0.1-0. 2 ray leigh. The adjustment to the linear factors F - and n, Y required 
to make the model predictions fit these data is about. .ten times larger 
than that required by any other combination of predictions and data, 
including the cell-empty case for the same interval. Therefore, we conclude 

O 

that these cell-full counts were not caused by 584 A flux. Instead, we 
believe the instrument observed several counts per second of background 

j 

flux of unknown, wavelength. We cannot determine the number of rayleighs 
of this radiation present, since the calibration of our detector varies 
with wavelength. - 

We then proceeded with the fitting of the data to the model predictions . 
First, we subtracted from each data point a constant 2.5 counts per second, 
to correct for the background count rate detected. Second, we converted 

the resulting count rates to rayleighs of flux by dividing by the instrument 

o -X -1 

584 A sensitivity of 18 counts s rayleigh . Finally, for each discrete 

interval of data collection, we determined the appropriate factor by 

which the nominal product F q must be multiplied to best fit (minimum 

chi-square) the envelope of predictions to the observed fluxes. We call 

this quantity the scale factor. The observed fluxes, with background 

removed, and the scaled model predictions, are displayed in Figures 1-2, 

which correspond to the eight distinct intervals of data collection which 

we have discussed. The envelopes of the model predictions are indicated by 
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shading, Solid lines, with representative error bars, depict the observed 
fluxes. The error bars do. not reflect calibration errors of the instrument, 
since these presumably result only in an unknown, constant correction to 
the nominal instrument sensitivity of 13 counts s ^ rayleigh which 
will automatically be factored out in determining the scale factors. 

The specific scale factor used for each Figure is displayed below 
the ordinate, and measures the amount by which the product of the solar 
584 A flux F q and the interstellar neutral helium density n^ assumed in 
the model differed from that actually present. Thus a scale factor of 
0.5 means that the product F^ actually present was half that assumed 

by the model. The scale factors vary from 0.63 to 0.84. 

The only background allowed for in our data reduction was a Poisson- 
distributed random term with constant mean rate 2.5 counts s If 
variations in the background sources mentioned caused that mean rate to 
change during the intervals of data collection, then the Poisson variances 
which we have used in calculating chi-square are too small. Variations 
in mean background rate of as little as 3 counts s ^ would increase the 
variances by factors of 1.5 to 4 and thereby reduce the calculated chi- 
squares by about a factor of two. The background sources just discussed 
are more than sufficient to account for such variations, although their 
quantitative behavior is insufficiently known to make more precise 
calculations of their effects on chi-square . 

In light of these background sources, we believe that our new model 
reasonably describes the HGD data. 
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The average scale factor for the eight distinct intervals of data 

collection is 0.76 ± 0.06..-. Thus, if our model predictions are to agree, 

on the average, with the reduced observations, we must adjust the values 

7 -5 -1 

of n^ and F^ so that their product is 0.76 x 10 photons cm s 

There has been -much discussion of the correct average value of the solar 

584 A flux F . Hall and Hinteregger (1970) reported F q = (8.9 ± 1.8) x 10 8 
-2 “1 

photons cm s and Hinteregger (T9 76) has recently revised this value 
upward by 30-60 per cent. Rusch et al. (1976) have tabulated Atmospheric 

g 

Explorer C observations for 1974 June 8, showing F q = (15.5 ±4.7) x 10 
-2 -1 

photons cm s . Maloy et al. (1975), observing on 1974 November 27, found 

• O ^-t 

F q = (26 ± 3) x 10 photons cm s , and have later (Carlson 1976) 

revised the uncertainty to 30 per cent. Heroins and Higgins (1977) 

8 ~~2 “X 

reported values of F q between 13 and 15 x 10 photons cm s , on five 

separate occasions from 1971 to 1976. To encompass these diverse results, 

8 2 —X 

we shall adopt F q = (20 ± 10) x 10 photons cm s as an appropriate 

average value for the remainder of this discussion. 

7 -2 -1 

Since the product n^ e F^ must equal 0.76 x 10 photons cm s , we 

—3 —3 

conclude that n He = (3.8 ± 2 .1) x 10 cm . The error given is the 

root sum square of the error assumed in F q , the calibration error of the 

instrument, and the standard deviations calculated in averaging the scale 

factors. Consequently, three major conclusions follow from the analysis 

° 

of the Apollo-Soyuz nighttime 584 A helium observations. First, tbe total 

. O 

solar 584 A flux may vary by as much as 20 per cent on a time scale as 
short as an hour or less. Second, the data and analysis generally support 
most recent conclusions about local ISM parameters derived from observations 
of resonantly scattered solar EUV flux. More specifically, our observations 
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are consistent with the interstellar gas moving toward a = 72°, fi = +15 
at a distant bulk speed of 20 km s \ with a temperature of 5000-20,000 K. 

Third, the neutral helium number density in the local ISM is (^3.8 ± 2.1) 

-3 -3 

X 10 cm . 

O 

It is clear that future interpreters of 584 A observations of the local 
ISM must be most cautious in drawing conclusions based on exact values of 

O 

the solar 584 A flux without near-simultaneous full-disc measurements of 
that quantity . 

The fact that our model fits were obtained with values of the ISM 
parameters generally similar to those recently used by other authors is 
more than a mere confirming. ^observation, for as we have discussed above, 

O 

most previous analysis of helium 584 A observations has been by means of 
a model strictly valid only for very low ISM temperatures. Our model does 
not possess this restriction, and so the broad agreement of our analysis 
with those which preceeded it reduces the possibility that the previously 
obtained fits of a low-temperature model to a high-temperature gas may have 
in some way been specious or misleading. 

The value of which we report is consistant with more recent 

-3 

conclusions. Freeman et al. (1977) concluded n„ < (0.004 ± 0.0022) cm 

on the basis of observations made during the Apollo-Soyuz mission, with a 

different instrument. Fahr, Lay and Wulf-Mathies (1977) have supported 

this conclusion. Given the nominal 0.085 relative cosmic abundance of helium 

atoms to hydrogen atoms (Allen 1973a)', our value for n He implies a -number 

-3. 

density n of neutral hydrogen in the nearby ISM of 0 .045 ± 0 .025 cm , 

H 
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2.2 The Daytime He I, 584 A Emission. 

Some of the MA--088 channel il 3 data were taken while the Sun was well 
above the horizon. These ’data have been examined to see whether the 
daytime count rates, which are much higher than the nighttime rates, can 
be explained by airglow lines of helium or of other atmospheric species. 

The results of this analysis are described in detail in a paper prepared 
for publication in the Journal of Geophysical Research and appended to this 
report as Paper 3. 

We imposed several additional constraints on the data selected for 
analysis and discussion for this case. We required that the line of sight 
point within 97° of the local zenith, which corresponds to the lower edge 
of the field of view being tangent to the Earth’s limb. We excluded data 
taken in or near the South Atlantic Anomaly. And to avoid contamination 
by sunlight reflected from the spacecraft, we included only data taken 
when the line of sight pointed more than 60° from the Sun. These 
conditions were met by a total of 1105 seconds of data, all taken on 
July 20 or 21, 1975, from an orbital altitude of about 225 km. 

With gas absorption cell empty, the count rates detected ranged from 
81 to 1880 counts s . With the cell filled, the corresponding range of 
rates was 12 to 811 counts s We correlated the reduced cell-full and 
cell-empty data with several parameters upon which airglow intensities 
might depend, including solar zenith angle, line-of-sigh't elevation, and 
angle between line of sight and the Sun. 

The only correlation we found was with solar zenith angle, which is 
also the angle in arc on the Earth’s surface between the shb solar point 
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and the spacecraft position. In figure 3 we show the cell-empty and cell- 
full count rates plotted against solar zenith angle, with typical error 
bars. The straight lines are linear regressions fits for the two cases. 

The r-square parameters for the two fits are 0.686 for the cell-empty data 
and 0.195 for the cell-full data. This dependence of the observed fluxes 
on spacecraft position is qualitatively consistent with the inodel of He X 

O , * 

584 A scattering from an optically thick neutral helium geocorona, developed 
by Meier and Weller (1972). 

As shown in Figure 3 , the cell-full count rate tracked the cell-empty 
count rate remarkably well. Although the fluxes observed varied by an 
overall- factor of twenty, the cell-full count rate remained at about 0.3 of 
the cell-empty count rate throughout. 

This behaviour of the HGD data has permitted us to come to the following 
conclusions : 

first, by subtracting cell— full from cell-empty count rates, we 
conclude that the daytime He I EUV geocoronal flux intensity spans the 
range 2-70 xayleighs. This value compares favorably with the 25 rayleighs 

Cv 

of 584 A flux reported by Christensen (1976) between 210 and 240 km at a solar 
zenith angle of 74°. We stress that because of the use of the gas 
absorption cell, our measurement is not contaminated by any radiation 
whose wavelength differs from the He I ground-state lines by more than a 

O 

few mA. We also restate that our data set is limited by several. constraints 
— notably, the line of sight always points 60° or more from the Sun. 

Second, inasmuch as the gas cell was designed to absorb all geocoronal 



- 12 - 


584 A flux, the substantial 500-7Q0 A flux observed through, the filled 

O 

cell must have soma other origin. Interplanetary 584 A flux might be 
Doppler-shifted sufficiently to be transmitted by the gas cell, but the 
theoretical work cited indicates that this flux is much too weak to explain 
our cell-full data. We cannot calculate an absolute intensity for this 
radiation, since its wavelength or wavelengths, and hence the correct 
calibration factors for our detector, are not known. 

Third, the strong correlation between the cell-full flux and such a 
local parameter as spacecraft position, suggests that this flux originated 
near the Earth, The intensity of this radiation varies with the observer’s 

O 

position in a manner similar, to that of the 584 A flux, at least for the 
range of observation, positions and directions considered here. 

Carruthers and Page ( 1972 , 1976) and Christensen (1976) have reported 
possible detections of dayglow spectral lines between 600 and 700 A 

O 

With intensities comparable to He I 584 A. Any of these lines might 
explain our cell-full fluxes. More detailed investigations with high- 
resolution spectrometers will be needed to define conclusively the spectral 
distribution of the radiation in this band. 

3. A HOT MODEL OF THE INTERSTELLAR MEDIU M 

In this work the distant ISM is described by a Maxwell-Bo It zmann 
distribution of temperature T, moving with bultc speed V- in a certain 
direction. This initial condition makes .the analytical description of the 
phase-space number density function n(v,s) a much more difficult problem. 
Furthermore, the use of such numerical techniques as Monte Carlo methods 
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to evaluate n(v,s), requires large amounts of time on even the fastest 
of present-day computers. It is therefore not surprising that most 
workers in the study of the local ISM continued to use the modified cold 
model, even when early results indicated that thermal speeds might in 
fact not he small. 

In 1975 Wallis pointed out that an analytic solution for n(v,s) under 
the initial conditions just given had been obtained by Danby and Camm 
in 1957, in the context of a different astrophysical problem. The Danby 
and Camm expression for n(v,s) contains an exponential of a complicated 
trigonometric function which is difficult to manipulate' analytically, 
furthermore, although Dirac delta functions in a(v , v, v.) allow two of 

* , i 

the integrations to be performed immediately, the calculation of an 
observed flux using the Danby and Camm n(v,s) nevertheless requires a 
four-dimensional Integration, three dimensions more than with the modified 
cold model. 

We have developed a program to evaluate the scattered intensity by 

numerical integration, using the Danby and Camm number density function, 

2 

modified by the inclusion of the factor exp (-B r_, Q/Vu) which was used by 

o i, 

•preceding authors (e.g. Blum, Pfleiderer and Wulf-Mathies 1975) to allow 
for the loss of gas atoms by ionization due to solar EUV. In this factor, 

8 is the rate at which gas atoms are ionized at one astronomical unit r 
from the Sitn, 0 is the angle through which a gas atom has orbited about 
the Sun, and p is the impact parameter of the orbit in question. We refer 
to this model as the hot model. 

Our hot model uses Simpson’s rule to perform all four integrations. 
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In testing the program, we varied the number and spacing of points at 
which the integrand was evaluated, until we observed an apparent asymptotic 
limit to the calculated value of the scattered flux I. For all further 
work, we used a number and spacing of points which produced a value of I 
within ten per cent of the observed limit. We checked the program validity 
by using several limiting cases singlyyand in combination, such as setting 
solar mass and radiation pressure equal to zero, letting AA approach 
infinity, and allowing T to trend toward zero. Calculations using the 
program also display symmetries reflecting the geometrical symmetries of 
the physical situation involved. 

The program provides the total observed flux and the observed line 
profile. It has provision for also calculating the flux and line profile 
observed when a helium gas absorption cell is interposed between the 
observer and the line of sight. We have so far used the program to 
investigate the behavior of neutral helium as observed via the resonant 
scattering of solar EUV at 584 A as it pertains to our HGD data. Each 
evaluation of the integral (1) requires . three to fifteen minutes of central 
processor time on a CDC 6400 computer. 

Fourteen physical parameters determine -the predictions of the models. 
These are three parameters describing the bulk velocity of the distant gas, 
the distant gas temperature T and neutral helium number density n , the 

* CO 

rate g at which the gas is ionized at one A.U. from the Sun, the FWHM AA 

of the solar resonance line, the solar line flux F at one A.U. from the 

3 o 

Sun, the ratio y of the pressure of light to the force of solar gravity 
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on a gas atom, three parameters characterizing the position of the 
observer and two describing the line of sight (LOS) of the observation. 

The results of this work demonstrate that when the modified cold 
model and the newly developed hot model predictions for helium atoms are 
evaluated for physical parameters in the range presently believed correct, 
the hot model predicts broader intensity maxima and smaller total 
yariations in flux . level tban does the modified cold model. The shapes 
qf the intensity maxima are sufficiently complicated that the difference 
in shape is difficult to quantify, but our calculations indicate that 

_ Q 

the patterns of 584 A glow predicted by the two models are markedly 
different within 60 to 90° away from the projection of the downwind 
symmetry axis on the celestial sphere. We conclude that the modified 
cold model should not be used to analyze observations of this region. 

At greater distances from the downwind axis, the shapes and spacings 
of the isophotes predicted by the two models are more similar, .although 
the hot model predicts several times as much flux as the modified cold 
model. The modified cold model may possibly be of use for interpreting 
observations of this part of the sky, provided the analysis used depends 
only on the relative intensities of flux observed from different points . 

O 

Future observers of 584 A helium glow should address the question of 
model validity more specifically. 

The hot model predicts substantial variations of observed flux with 
temperature, on lines of sight for which the modified cold model predicts 
no variation at all. Furthermore, even when the modified cold model does 
predict a variation of flux with temperature, the hot model predicts 
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variations of a different kind. The use of the' modified cold model to 
determine the precise temperature of the local ISM therefore appears 
invalid, 

O 

The conclusions drawn by previous observers of 584 A flux who have 
used the modified cold model to interpret their data must be reassessed. 

The hot model does not fit the data of Weller and Meier (1974) as well 
as the modified cold model did. We therefore suggest that further study 
pf both model and data is in order, to determine the cause of the discrepancy 
and the actual values of the ISM parameters that the data imply. When the 
hot model is used to analyze the data of Freeman et al. (1976) the derived 
lower limit to gas bulk speed is slightly rearer to the values of 
obtained by other means than was the case with the modified cold model. 

The theoretical question of the behavior of the local ISM is far from 
closed. Both the modified cold model and the hot model derive from two 
assumptions. First, the flow of the ISM through the solar system is 
assumed to be collisionless, so that each atom proceeds along a hyperbolic 
orbit determined by celestial mechanics. Second, the nearby ISM is assumed 
sufficiently tenuous that only singly-scattered photons will reach the 
observer. Wallis Q.973, 1974, 1975) has questioned both of these assumptions. 
Further analytic and numerical study of their validity may be necessary before 

O ' 

experimenters can have complete confidence in the interpretation of 584 A 


data. 
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FIGURE CAPTIONS 


Figures 1- 


Figure 9 : 


Reduced observed fluxes and model predictions for the 
eight discrete intervals of data collection. Open 
circles = cell-empty fluxes. Solid circles = cell-full 
fluxes. Representative error bars are shown. Shaded 
area = envelope of model predictions as temperature varies 
from 5000 to 20,000 K. The model predictions have been 
adjusted so that the product F q - A x 10 ^ photons 
cm ** s where the value of A is displayed at the bottom 
of each Figure. • 

Cell-full (lower) and cell-empty (upper) count rates 
versus solar zenith angle. Calculated linear regressions 
fits shown by diagonal lines. 
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ABSTRACT 


A photometer sensitive at the 584 line of HeT, incorporating a helium 
gas resonant absorption cell, was flown on the Apollo-Soyuz Test Project 
in July 1975. The instrument observed much of the nighttime sky, and 
returned 42 minutest" useable data. The data were analyzed by fitting to 

O 

a model of resonant scattering of solar 584 A flux from nearby interstellar 

helium, which is valid even when the interstellar gas is hot. Good model 

fits were obtained for an interstellar gas bulk velocity vector pointing 

-1 

toward os = 72°, 6 = +15°, with speed 20 km s with interstellar medium 
temperatures from 5000 to 20,000 K and with neutral interstellar helium 
density 0.0038 + 0.0021 cm~^. The results are consistent with the Sun 
lying in the warm, partially ionized periphery of a cold interstellar cloud 
surrounded by old supernova remnants, as described by McKee and Ostriker 


(1977). 



1. INTRODUCTION 


Observations of interstellar gas passing .through the solar system 
provide information about the temperature, number density, composition 
and bulk velocity of the local interstellar medium (ISM). . A number of 
authors (Blum and Fahr 1970; Holzer and Axford 197.1; Axford 1972; Johnson 
1972; Paresce and Bowyer 1973; Wallis 1973, 1974, 1975; Meier 1977; 

Freeman and Paresce 1978) have described how neut^cal atoms passing near the 
' Sun are gravitationally focused into a densf ty distribution whose shape 
depends. in a complicated manner both on physical parameters of the ISM and 
on the local effects of photoionization b} T the Sun, collisional ionization 
and heating by the solar wind, collisions with other atoms, and the pressure 
of sunlight. Study of the details of this distribution allows .untangling 
of the various parameters by model fitting. Studies of the distribution 
of the local ISM have been made by observation of solar ultraviolet and 
extreme ultraviolet (EUV) light scattered by hydrogen (Bertaux and Blamont 
1971; Thomas and Krassa 1971; Thomas 1972; Bertaux et al. 1976), and helium 
(Paresce, Bowyer and Kumar 1973a, b, 1974; Weller and Meier 1974; Freeman 



et al. 1976a, .1977; Ajello 1977). These experiments have established 
the direction toward which the local ISM moves' as within about fifteen 


■degrees of a = 72°, 6 = +15°, and have determined the bulk speed of the 

distant gas to be within a few km s”-* - of 20, but there remains a somewhat 

greater controversy about the temperature and number density of the gas. 

There has been a substantial theoretical difficulty with the models 

used to interpret the data from these experiments. Most observers of helium 

have performed their model-fitting using a model which is strict3.y 

valid only if the thermal speeds of individual atoms are very small compared 

to the bulk speed of the distant gas. At a temperature of 10^ K, roughly 

the mid-point of the range discussed in the references cited, the actual 

helium thermal speed is 4.6 km whereas the bulk speed V seems to be 

is 

about 20 km s~^, as discussed above. There is therefore room for doubt 


about the applicability of the low-temperature models and the accuracy of 
conclusions derived by using them. Models- valid at higher temperatures are 
analytically complicated and require large amounts of computer time to 


evaluate, and have only recently become available (Meier 1977; Freeman and 



Paresce 1978) . In the present work, we use the model of Freeman and Paresce 


O 

to interpret data obtained with a 584 A photometer, equipped with a helium 
gas resonant absorption cell, during the 1975 Apollo-rSoyuz Test Project. 


II. INSTRUMENTATION 


Complete details of the construction and calibration of the Apollo- 

O 

Soyuz helium 584 A photometer have been published elsewhere (Bowyer et al. 
1977). In brief, the instrument consisted of a channel electron multiplier 
(CEM) whose optical axis traversed a gas cell 10 cm long. The cell was 
sealed by metal foil EUV filters, one tin and one aluminum, each approximately 

O 

1200 A thick. The cell was alternately filled with helium at a pressure 

of 1.4 torr for a period of approximately 2.5 s, and then vented to space 

for an equal interval. ThjjCelium column density within the filled cell 
17 -2 

was 4.9 x 10 atoms cm . The helium column absorbed more than 99 per cent 

o 'o 

over a band about 40 mA full width, centered on the 584 A line. This 

absorption feature was sufficiently broad to absorb solar 584 A flux 
\ 

-scattered by neutral helium in the Earth's upper atmosphere, even at several 
times the nominal (Jacchia 1971) exospheric temperature of 10^ K, and even 



when that scattered flux was Doppl'er-shifted by the entire 8 km s''*' 
orbital speed of the spacecraft. 

When the cell was vented, the partial pressure of helium within it V7as 

calculated to fall to less than 10“^ torr in a few tenths of a second, which 

11 -2 

corresponds to a helium column density not exceeding 3 x 10 atoms cm' . 

At this column density the absorption feature would have had a full width 

O 

at half maximum of 2 mA and a maximum absorption, at line center, of 12 per 
cent of the incident flux. Such an absorption feature is narrow compared 

O 

with the calculated profiles of the solar 584 A line scattered from the 
-local ISM. Hence absorption caused by the vented gas cell was negligible. 

An end-to-end photometric calibration of the instrument was performed 
before flight, using emission lines of hydrogen, helium, neon, and argon at 

O 

wavelengths of '200 to 2000 A. NBS-calibrated A^O^ photodiodes and gas 
proportional counters were used as reference standards. With gas cell vented, 

, o _1 

the instrument sensitivity to diffuse 584 A flux was 18 + 3.6 counts s 
rayleigh - ^. The product of CEM quantum efficiency and metal foil filter 
transmissions was 4 x 10 , 


and the field of view was approximately 14 degrees 



in diameter. In the laboratory, the instrument exhibited a dark count rate 
of less than 0.1 xe count s A comparable count rate was obtained in 
flight with the instrument door closed. 

III. DATA 

The Apollo spacecraft carrying our instrument was launched in July, 

1975, into a near-circular orbit at 220 km altitude. The orbital plane 

was inclined at 55° to the celestial equator, and approximately contained the 

Earth-Sun line. Data were taken at various times during the mission while 

the spacecraft was maneuvered so that the instrument field of view surveyed 

most of the celestial sphere. The instrument was also turned on at other 

times to acquire additional data from random lines of sight. The instrument 

was not operated during or immediately after purges of gas or liquid from 

the spacecraft, ox when the reaction control system rocket motors nearest 

the instrument were activated. Approximately ten hours of data were obtained. 

* / 

NASA provided detailed spacecraft aspect and trajectory information for this 
time period. 

Inspection of the data revealed that the observed count rates were 


100-1000 times as large when the Sun was above the horizon as when it was 



below the horizon. We ascribe these high daytime count rates to scattering 
of solar EUV radiation from residual atmosphere at and above the spacecraft 
orbit, and have discussed the matter in detail elsewhere (Freeman et al. 1978). 
To make certain that our conclusions about flux scattered from the local 
ISM were not affected by this daytime radiation, we restricted the present 
analysis to data taken when the Sun was at least 30° below the-horizon. Since 

the Earth's limb is depressed 15° below the horizontal when viewed from 220 km 

* 

altitude, this constraint corresponded to a solar zenith angle of at least 
135° . To minimize the affect on our analysis of errors in calculating 

O 

absorption of 584 A flux by the atmosphere above 220 km, we chose to analyze 
only those data for which the zenith angle of the line of sight (LOS) was 
less than 60° . • Finally, we excluded .from analysis data taken when the 
spacecraft was in the South Atlantic Anomaly. After these three constraints 
were imposed, 53 minutes of data remained for study. 

While these d^ta were being collected,^ the LOS scanned the sky at a rate 
of at most 0.5° s“^. We therefore deemed it appropriate to bin the data 
into 30-s intervals, corresponding to the time necessary for the LOS to 



traverse at most one field width. During each. such interval we separately 


treated data taken when the gas absorption cell was full and empty, and 

ignored data taken while the cell pressure was changing. In order to 
' ' \ 
eliminate spikes of extraneous noise and telemetry dropouts from the data, 

we subjected each bin of 0.1-s data samples to a filtering process in which 

the mean and standard deviation of each set of samples was calculated, and 

then all samples more than three standard deviations away from the mean 

were discarded. The process was repeated until no more samples were discarded 

The amount of data threrby eliminated was less than five per cent and was 

consistent with Poisson statistics. The resulting count rates for the various 

bins ranged between 2 and 60 counts s - ^. The standard deviations obtained 

from the statistics of the individual 0.1-s data samples were consistent 

with those expected from Poisson statistics, and were generally in the range 

— 1 

0.5 - 1.5 counts s . 


We corrected the observed count rates for atmospheric absorption of 

O 

584 A flux along the lines. of sight by using the atmospheric_„models of 
Jacchia (1971). The solar activity parameters required as inputs to, the 



models were taken from Lincoln (1975). We used the 584 1 absorption cross- 

sections of Hall, Schweizer, and Hinteregger (1965) for N , 0 and 0 . Our 

2 2 

calculations were based on the spacecraft position and LOS orientation at 
the mid-point of each bin. The effect of these corrections was to increase 
the observed flux levels by 10 to 30 per cent. We have not included in our 

o 

stated errors the effect of any errors in the Jacchia models or the 584 A 
absorption cross-sections of K 2 , 0 and O 2 . 

O 

It was also necessary to allow for resonant absorption of 584 A flux by 

■t •/ r\ 

■the 1-2 x 10 x atoms cm~ (Jacchia 1971) of atmospheric helium remaining 
above the spacecraft. This correction is model-rdependent , for it depends 

O 

on the spectral profile of the 584 A line scattered from the ISM helium. 
Therefore, we chose to incorporate this correction into the models used for 
predicting the scattered lines, which we shall describe in the next section. 

We investigated the .possibility of other terrestrial phenomena affecting 
the data by correlating the binned, absorption-corrected data with LOS zenith 
angle, solar zenith angle, LOS-Sun angle, local strength of the geomagnetic 


field, magnetic latitude, and component of spacecraft velocity along the LOS. 



The only correlation we. found was with magnetic latitude. In Figure 1 we 

/ 

show cell-empty count rates plotted against magnetic latitude. There is a 
pronounced increase in count rates south of magnetic latitude -r5Q°, and a 
noticeable peak within 10° of the magnetic equator. These increased count 
rates may be due to localized airglow, particle precipitation, or some 
other cause. We decided to restrict the data under analysis further, by 
excluding data taken south of magnetic latitude -50° or within 10° of the 
magnetic equator. After this additional exclusion there remained 42 minutes 
of data. In Figure 2 we depict by hatching the portions of the sky viewed 
while these data were gathered. 

We encountered one additional uncertainty in reducing the data. As 

O 

discussed by Schmid tke (1976), the total solar 584 A flux varies from day to 
day by as much as thirty per cent. This short-term variability appears to be 
sufficiently unsystematic that we ar*e reluctant .to interpolate between Schmidthe’ 
published data points, which are spaced at least 24 hours apart. We are not 

t 

O 

aware of any studies of solar 584 A variability on shorter time scales. 

— — The 42 minutes of Apollo-Soyuz data we have selected for analysis are 



composed of eight separate intervals, 1—10 miliutes long, spaced 1 to 24 hours 

* o 

apart. To allow for possible solar 584 A flux variations from interval to 
interval, we decided to allow the assumed value of the flux to vary, from 
interval to interval, while searching for model fits-. 


IV. MODELING OF PHOTOMETRY OF THE LOCAL ISM 


A general formalism describing radiation resonantly scattered from 
hydrogen and helium in the local ISM, in the approximation that thermal 
speeds of individual atoms are small compared to the bulk speed of the gas 
as a whole, is well established (e.g. Freeman et al. 1976a and references 
therein). We have recently described (Freeman and Paresce 1978) a revised 
model which does not require the approximation of small thermal speeds. The 
model evaluates the following integral for the observed flux I: 

3 : = CeGO;) f Gi>, , .s') ^ (-0*. , ^ , "0;^ p CoC-sS)/! Gj - ,G) 

c{ V/ 

In this integral, X is the number of photons detected per second by the 

instrument, v. is the radiation frequency in an inertial frame attached to 

A- 

the Sun, s is path length along the LOS, and v is gas velocity. The quantity 
E(v^) is instrument sensitivity as a function of wavelength, F(v A ,s) is the 
solar flux at the point described by s, a(v , v, v ) is the total cross- 



section for an atom moving with velocity v to scatter light of frequency 
v A at the Sun in such a manner that the observer perceives its frequency 
as 0(s) is the angle at s between the Sun and the observer, p (0(s)) 
is a phase function, and n(v^s) i's the gas phase-space number density at s. 
Our revised model differs from those which preceded it (except for the 

A 

1977 model of Meier) in that we use for n(v,s) an analytic expression due 

to Danby and Camm (1957), which is valid even for high temperatures of the 

distant ISM. Although Dirac delta functions in a(v A , v, v^) allow two of 

the integrations in equation (1) to be performed immediately, numerical 

solution of the remaining four-dimensional integral requires 3 to 15 

minutes of central processor time on a CDC 6400 computer. Our program 

simultaneously evaluates equation (1) for two different functions E(vi'), 

corresponding to gas absorption cell full and empty . 

✓ 

The predictions of the model depend on fifteen physical parameters, in 

addition to the function E(v.). The parameters describe the Bulk speed V , 

x D 

temperature T, helium spatial number -density n and direction of motion of 

He 


the distant gas, the position and LOS of the observer, the component of the 



observer’s velocity along the LOS, the total flux F of the solar 584 A line 

o 

as measured at 1 AU from the Sun, the FWHM AX of that line, the rate g 

' o 

at which neutral helium at 1 A.U. from the Sun is ionized by solar 

radiation at wavelengths less than 504 A, and the ratio y between the 

pressure of sunlight and the force of solar gravity on a neutral helium atom. 

For the modeling used in the present work, we assumed that the local 

ISM flows toward the direction a = 72°, 6 - 4-15°, given by Weller and Meier 

(1974). Previous observations of the local ISM with gas absorption cells 

have placed limits on the gas bulk speed V . Bertaux et al. (1976), observing 

B 

hydrogen, found between 12.5 and 21 km s“‘*‘, and most likely between 18 and 
21 km s“ , Freeman et al. (1976a), observing helium, found a lower limit of 
10-15 km s'* 1 . Adams and Frisch (1977) have used Copernicus higher e solution 
spectrometer scans of Lyman-alpha background to measure the speed by Doppler 
.shift, and found about 22 km $~ 1 . We therefore- chose = 20 km s _1 for our 
modeling. 

» 

Weller and Meier (1974) and Meier (1977) used AX = 0.12 A and B o = 6.8 
—8 \ 

x 10 s Delaboudiniere and Crifo (1975) summarize recent work indicating 

o 

that AX is between 0.11 and 0.15 A. Banks and Kockarts (1973) indicate that 



g = (8.0 + 4.0) x 10"~^ s^ 1 . We selected AA == 0,12 A. and $ « 6,8 x 10"^ s ^ 

o — o 

for the present work. We chose y = 0, as indicated by Fahr (1974), and 

we used values for the observer’s position, velocity and LOS appropriate 

to the mid-points of the various data bins. 

There is much uncertainty in the correct value of n^, and ' as we have 

indicated, F q is known to vary substantially on a short time-scale. However, 

both of these quantities appear in the evaluation of (1) only as overall 

linear factors. Therefore, it is not necessary to justify the values used 

in. the initial modeling, because the model predictions need only be scaled 

appropriately at a later stage in the analysis. We selected the numerically 

o 9 —2 —1 

convenient values = 0.01 cm“ J and F q = 10 cm s for inputs to our 
computer program. 

The work of Paresce, Bowyer, and Kumar (19/4), Weller and Meier (1974), 

and Bertaux et al, (1976) suggests that T is within a factor of two or so of 

» 

10^ K, Therefore, we evaluated our model for values of 5000, 10,000 and 
20,000 K, • 

O 

Finally, we corrected the 584 A line “profiles predicted ; by the model 



for resonant absorption by atmospheric helium above the spacecraft, using 


the atmospheric models of Jacchia (1971) . The' effect of this helium was 
to impose on the scattered lines an absorption feature typically 12 mA 
FVJRM, which absorbed more than 99 per cent of incident flux over a typical 
bandwidth of 10 mA. The feature was Doppler-shif t-ed by the relative motion 
of the spacecraft to the atmosphere, but it was always contained within the 
absorption baiid produced by the filled helium gas absorption cell, and 

therefore could not alter the model predictions for - cell-full flux, 

, * 

The moddl predictions for the data discussed here were relatively 

insensitive to temperature, typically varying by less than 10 per cent as T 

was raised from 5000 to 20,000 K« This property is a straightforward 

consequence of the geometry of our observations. The lines of sight passed 

through only regions well upwind of the Sun, where the effects of temperature 

are known to be small (Freeman and Paresce 1978).. 


V. COMPARISON OF I)ATA AND MODEL PREDICTIONS 

We evaluated the model predictions for the physical parameters 
corresponding to the mid-points of the various data bins. In all cases. 





the cell-full count rates and predictions vex;e. substantially less than 
the corresponding cell-empty values,' Almost all cell-empty flux predictions 
increased with temperature, as did most cell-full predictions. 

t thaS<L 

We studied the non-584 A "background count rate by inspecting^observations 

O 

for which our 584 A models predicted a count rate of less than a few tenths 
of a count per second. During about 20 s of such observations, the background 
rate detected was 2.5 counts 

We then proceeded witk the fitting of the data to the model predictions, 
First, we subtracted 2.5 counts per second from each data point to correct 
for the background rate. Second, we converted the resulting count rates to 

O 

rayleighs of flux by dividing by the instrument 584 A sensitivity of 18 
counts s'* 1 rayleigh"'-*- . Finally, for each discrete interval of data collection* 
we determined the appropriate factor by which the- nominal product n^ e F q must 
be multiplied to best fit (minimum chi-square) the envelope of predictions 
to the observed fltixes, We call this quantity the scale factor. The observed 
fluxes, with background removed, and the scaled model predictions, are 
displayed in Figures 3 and 4, which display the eight distinct intervals of 



data collection which- we have discussed. The envelopes of the model 
predictions are indicated by shading. . Circ lei , with, representative 
error bars, depict th^bserved fluxes. The error bars do not reflect 
calibration errors of the instrument, since these presumably result only 
in an unknown, constant correction to the nominal instrument sensitivity 
of 18 counts s~i rayleigh - '-*-, which will automatically be factored out 
in determining the scale factors. We shall include the effect of sac,h 
calibration errors where appropriate in further analysis. 

The specific scale factor used for each interval is displayed below 
the ordinate, and measures the amount by which the product of the solar 

O 

584 A flux F and the interstellar neutral helium density n, assumed in the 
o He 

model differed from that actually present. Thus a scale factor of 0,5 means 

that the product n F actually present was half that assumed by the model, 

He 0 

The scale factors vary from 0,63 to 0,84, The scale factor used for each 
Figure is the average of the three separate scale factors obtained for 
best fits of the data to the model predictions for the three temperatures. - 
used. Each separate scale factor was derived by minimizing the chi-square 


statistic for the fit in question. We have scaled the predictions in the 


Figures using the average scale factor, rather than using for each prediction 


the specific scale factor which minimizes its own chi-square, in order 
. to display the variations in predictions with temperature more clearly, 
The data displayed in Figures 3b and 4b are .well -fit by the models. 

In Figure 3b, the minimum chi-squares for 5000, 10,000 and 20,000 K are 
respectively 2.69, 2,85 and 3,10, which imply, confidence levels (with five 
degrees of freedom) of 0,75, 0,72, and 0,68 for the three temperatures 
respectively. In Figure 4b, there are 20 degrees of freedom and the minimum 
chi-squares are 30.59, 23.12, and 10.20 respectively. The confidence levels 
in this case are 0,06, 0,’28, and 0,96 for 5000, 10,000 and 20,000 K 
respectively. 

In- all other figures, the chi-square fits are worse. Chi-square 
per degree of freedom ranges from 2 to 7, indicating very low confidence 

O 

levels for the models. Spatial variations in the non-584 A -flux may have 
been responsible for much of this discrepancy. For example, a 3 counts s“^ 

t 

variation in this flux would decrease chi-square by a factor of two. 

Possible sources of background counts include nighttime EUV airglow 

O 

lines (Carruthers and Page 1972), geocoronal helium 584 A radiation (Freeman 


et al. 1976), and the unknown source which caused .our high count rates 
near magnetic latitudes 0°-and -50% The spatial variations observed in 
these sources are more than sufficient to change the instrument background 
rate by several counts s“^« Unfortunately, the" quantitative details of 
these spatial variations are : insufficiently known to perm.it more precise 
calculations of their effects on chi-square* 

In light of these background ■ sources, we believe that our models 
reasonably describe the data presented here* With the reductions in chi- 
square just discussed, the data do not rule out any of the three temperatures 
we used as model parameters . 

The average scale factor for the eight distinct intervals of 
data collection is 0*76 + 0.06* Thus, if our model predictions are to 
agree, on the average, with the reduced observations, we must adjust the 

" 7 c -i 

values of n Re and F q so that their product is 0,76 x 10 photons cm s % 
There has been much discussion of the correct average value of the 
solar 584 A flux F^* Hall and Hinteregger (1970) reported F Q = (8,9 jh 1,8) 

o 9 

x 10° photons cnT" s~ and Hinteregger (1976) has recently revised this value 



upward by 30-60 per cent, Rusch et al, (1976) have tabulated Atmospheric 

6 

Explorer G observations for 1974 June 8, showing F q = (15,5 + 4,7) x 10 
—1 

photons cm s Haloy et al* (1975), observing on 1974 November 27, found 

8 1 

E -= (26 + 3) x 10 photons cm*" s" 1 , and have later (Carlson 1976) revised 
o 

the uncertainty to 30 per cent, Heroux and Higgins (1977) reported values 

3 2 —i 

of F Q between 13 and 15 x 10 photons cm z s , on five separate occasions 

i 

from 1971 to 1976, To encompass these diverse results, we shall adopt F q = 

o 7 1 

(20 + 10) x 10 photons cm _/ s~ as an appropriate average value for the 
remainder of this discussion, 

7 — S —1 " 

Since the product n F must equal 0,76 x 10 photons cm J s , we 

He o 

o 

conclude that n = (3.8 + 2.1) x 10 ^ cm , The error given is the root 
He — 

sum square of the error assumed in F , the calibration error of the 

o 

, ^ 

instrument, and the standard deviation calculated in averag ing the scale 

factors. 

VI. DISCUSSION AND CONCLUSIONS 

O 

Two major conclusions follow from the analysis of Apollo-Soyuz 584 A 


helium observations presented here. First, the data and analysis generally 



support most recent conclusions about local ISM parameters derived from 

observations of resonantly scattered solar EUYlflux, Second, the neutral helium 

• -3 t3 

number density in the local ISM is (3.8 +_ 2,1) x 10 cm , Furthermore, 

O 

it is clear that future interpretations of 584 A observations of the local 
ISM must be cautious In drawing conclusions based' on exact values of the 

O 

solar 584 A flux without near-simultaneous full-disc measurements of that 
quantity. 

The fact that our model fits were obtained with values of the ISM 
parameters generally similar to those recently used by other authors is 

O 

significant in that most previous : analysis of helium 584 A observations 
has been by means of a model strictly valid only for very low ISM temperatures, 
Our model does not possess this restriction, so the broad agreement of 
our analysis with those which preceeded it reduces the possibility that the 
previously obtained fits of a low-temperature model to a high-temperature gas 
may have in some way been specious or misleading. 

The low value of the number density of neutral helium in the local ISM, 

0.0038 Hh 0,0021 cm , is in conflict with early results derived from some 

© 

of the first observations of ISM 584 A flux. Paresce, Bowyer, and Kumar 



(l973a, b, 1974) gave n - 0,032 cm 3 , and Weller and Meier (1974) stated 

that n Re was within the range 0.009-0.024 cnf 3 , These early results depended 
the assumed value of F q (as ours does) , and the authors in question 

used values of that quantity approximately half the one we have chosen. 

Their value would be reduced by a factor of two if they had used our 

value for F . 

o 

The value of n which -v?e report is consistent with more recent 
He 

conclusions. Freeman et al. (1977) concluded n^ < (0,004 + 0,0022) erf 3 
on the basis of observations made during the' Apollo-Soyuz mission, with a 
different instrument. Fahr, Lay and Wulf-Mathies (1977) have reported data 
vrhich support this conclusion. Given the nominal 0,085 relative cosmic 
abundance of helium atoms to hydrogen atoms (Allen 1973a) , our value for 
n He iraplies a nuraber density n R of neutral hydrogen in the nearby ISM of 


0.045 + 0.025 cm -3 . 


A value of n of 0.045 is a bit low to be consistent v with the nominal 

H , 

-3 . o 

0.1 cm obtained by 1216 A observations of hydrogen in the local ISM 
(Fahr 1974). But as Fahr has pointed out, the interpretation of these 
observations is very model dependent, and there is considerable scatter in 


the reported values of n , 

H 



value of n of 0*045 is consistent with, hydrogen number densities 
H 

determined from column density measurements within a few tens of parsecs. 

Table X lists densities obtained from 21 stars within 100 pc of the Sun, 

including 16 values consistent with column-averaged n at or below* 0,05 cm 

, H 

In particular, Evans, Jordan and Wilson 0-975) observed a CMi, at a distance 
of only 3.5 pc, and found n^ between 0,015 and 0,03 cm*' 3 . Perhaps, as 
Freeman et al. (1977) suggested, this region of very tenuous gas includes 
the solar system. 

The preceding discussion suggests some possible conclusions about the 


structure and nature of the ISM near the Sun, McKee and Ostri her (1977) 

Vh# 


5 g 

have recently developed a multiphase ISM model which explains the 10 -10 K 
ISM temperatures inferred from interstellar 0 VI absorption lines (York 1974, 
Jenkins and Meloy 1974) by hot, tenuous gas in supernova remnants (alse \ 

v t * 

see Cox and Smith 1974), but which also allows for the presence of substantial 


volumes of partially ionized gas at near 8000 K, with neutral 

i 

hydrogen number density near 0,09 cm~ 3 . This gas is produced bleating and 


ionization of small, cold interstellar clouds which are embedded in the hot 



supernova remnants. Our results could be explained if the solar system 
lay within such a volume of gas. If so, then somewhere within a few pc 
of the Sun, one may find the cloud’s cold core, 'with T = 50-100 K, 
n^ - 10-100 cm~3, and radius 0,5-10 pc (McKee and Ostriker 1977), In 
view of the 'small number of measured hydrogen column densities to nearby 
stars, it seems possible that a small cold cloud could have escaped detection 
even so close to the Sun as a few pc. 

The clouds in the McKee-Ostriker theory are immersed in a network of old 
■ super nc^a. remnants, some of which may have been reheated by more recent 
supernovae. Supernovae induce substantial long-lived bulk motions of 
interstellar material radially outward from their positions (Chevalier 1977. 
and references therein) , On allowing for the motion of the Sun toward 
a - 271°, 6 = +30°, at 19.7 km s"" 1 (Allen 1973b), we find that with respect 
to the. local group of stars, the ISM approaches from a = 192°, 6 = -67°, at 
16 km s -1 , with uncertanties of aobut 15° and 3 km s'" 1 . These coordinates 
define a position within the Scorpius-Centarus association of early main- 
sequence' stars, an association which might reasonably be expected to have 



produced several supernovas earlier in its lifetime^ The strong 0 VX 
absorption (log ^ = 12,96 cm ) reported by Jenkins and Meloy (.1974) 

to g Cen A, at a distance of only 81 pc, may support the hypothesis of a 
nearby old supernova remnant. 

Whatever the result of these conjectures, it is clear that observations 
of interstellar gas passing through the solar system can provide valuable 
insight into the structure and properties of the interstellar medium out to 
tens of parsecs from the Sun, 
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TABLE I 


Star 

Distance (pc) 

H 

."N 

O 

3 , 

! 

Reference for 
Column Density of H-I 

a Cen 

1.34 

0.20 ± 0.05 

Dupree et al. (1977) 

e Eri 

3.3 

< 0.12 

McClintock et al. (1976) 

a CMi 

3.5 

0.015 - 0.030 

Evans et al. (1975) 

3 Gem 

10.7 

<0.15 

McClintock et al. (1975) 

a Boo 

11.1 

< 0.1 

Moos et al. (1974) 

a Aur 

14 

0.03 ± 0.01 

Dupree et al. (1977) 

a Tau 

20.8 

< 0.2 

McClintock et al. (1975) 

a Leo 

22 

0.02 

Rogerson et al. (1973) 

a Eri 

28 

0.07 

Rogers on et al. (1973) 

a Gru 

29 

0.09 - 0.18 

York (1976) 

n UMa 

42 

0.005 

York (1976) 

cr Sgr 

57 

< 0.17 

Bohlin et al. (1977) 

HZ 43 

62.5 

< 0.013 

Auer et al. (1977) 

a Pav 

63 

< 0.01 

Bohlin et al. (1975) 

a Sgr 

80 

< 0.12 

Bohlin et al. (1975) 

3 Cen 

• 81 

0.13 

Bohlin et al. - (1977) 

3 Lib 

83 

0.06 - 0.13 

York (1976) 

£ Cen 

83 

< 0.39 ' 

Bohlin et al. (1977)' 

a Vir 

87 

0.037 

Bohlin et al. (1977) 

Feige 24 

90 

< 0.011 

Margon et al. (1976) 

X Sco 

100 

< 0.078 

Bohlin et al. (1977) 




FIGURE CAPTIONS 


Figure 1: Observed cell-empty count rate, corrected for atmospheric 

absorption, as a function of magnetic latitude. Representative 
error bars are shown. 

Figure 2 : The' celestial sphere, indicatni<^ by hatching the portions 

O 

surveyed by the Apollo-Soyuz helium 584 A photometer, while 
gathering the 42 minutes of data analyzed herein. 

Figures 3 & 4: Reduced observed fluxes and model predictions for the 

eight discrete intervals of data collection. Open circles 

= cell-empty fluxes; Solid circles = cell-full fluxes. 

Representative error bars are shown. Shaded area - envelope 

of model predictions as temperature varies from 5000 to 20,000 

K, The model predictions have been adjusted so that the 

-7 -5 -1 

product n F = A x 10 photons cm s , where the value 
of A is displayed at the bottom of each panel. All fluxes 
are in rayleighs, and each data point represents one 30-s 
bin. 
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EXTREME ULTRAVIOLET DAYGLOW OBSERVATIONS 
WITH A 

HELIUM GAS ABSORPTION CELL 


Jay Freeman, Francesco Paresce, 
Stuart Bov/yer, and Michael Lamp ton 



ABSTRACT 


During the Apollo-Soyuz Test Project, an extreme-ultraviolet photo- 
meter with a cyclically operated -helium gas absorption cell oBserved the 
daytime sky from an orbital altitude of 225 km. When the line of sight 
pointed more than 60° from the sun, the instrument detected 2 to 70 
rayleighs of flux scattered from neutral geocoronai .helium at wavelengths 
from 504 to 584 £. The instrument also detected other radiation in the 
band 500-700 £ of- similar spatial distribution to the helium flux, which 
was definitely not due to the He I 584 £ spectral line. Possible sources 
of this radiation are discussed. 

INTRODUCTION 


The extreme ultraviolet (EUV) dayglow between 100 and 1000 £ has been 
observed by a number of authors (Donahue and Kumer 1971; Meier and Weller 
1972, 1974; Carruthers and Page 1972, 1976; Kumar, Bowyer and Lampton 1973; 
Carlson and Judge 1973). Fluxes of several rayleighs to several thousand 
rayleighs have been detected and ascribed to specific spectral lines. How- 
ever, poor EUV spectral resolution, in part due to the relatively broad band 
passes of available EUV filters, has made wavelength identification difficult 
and may possibly have prevented the isolation of closely-spaced lines. A 
spectrometer measurement (Christensen 1976) provided 20 £ resolution between 
550 and 1250 £, which is still inadequate to separate many possible dayglow 


lines. 
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We have recently flown a photometer sensitive at 500-700 X, in which 
a helium gas resonance absorption cell permits separation of the 584 R 
He I flux from other radiation in this band. The observation allows 
confirmation of the coarse-resolution spectrometer measurements of 584 & 
flux just discussed, and permits study of the spatial distribution of both 
the 584 R flux and of other, unresolved radiation in the band. 


INSTRUMENTATION 


The Apollo-Soyuz helium glow detector has been described in detail by 
Bowyer et al. (1977). In brief, it contained four EUV photometers, each 
with a 13 mm diameter active area and a 14° FWHM circular field of view. 

One photometer consisted of a channel electron multiplier with a helium gas 
resonant absorption cell which was equipped with one tin and one aluminum 
thin metallic film. The EUV bandpass of the instrument was defined by these 
films, which formed the entrance and exit aperatures of the absorption 

cell and transmitted from about 500 A to about 700 A. This bandpass 

- - ’ 2 
contained all neutral helium ground-state resonance lines, from the Is -ls2p 

transition at 584 & to the continuum edge at 504 R. The cell was alternately 

evacuated and filled with helium at a pressure of 1.4 torr, with a helium 

17 —2 

column density of 4.4 x 10 atoms cm . The filled cell superimposed on 
the 584 R line an absorption feature of 100 mR FWHM, which absorbed more than 
99 percent of incident radiation over a band 40 mR wide. A narrower absorp- 
tion featurewas superimposed on each shorter wavelength ground-state helium 
line. Every 0.1 s the number of counts detected was telemetered, together 
with the helium cell pressure and other information. 
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The instrument was calibrated in the laboratory before flight using 
a hollow-cathode discharge lamp operated with various gases to yield 
radiation in the 200 to 2000 8 region. An NBS-calibrated A^O^ 
windowless photodiode was used as a reference standard. The 
measured photometric efficiency of the instrument, versus wavelength, is 
shorn in Figure 1. The overall sensitivity of the detector to diffuse 
584 X flux was 18 ±3.6 counts s ^ rayleigh The detector response was 
constant as long as its photosensitive surface was directly illuminated by 
the calibration beam. When the beam reached the detector after 
one or more reflections from surfaces adjacent to the optical path, the 
response was reduced by a factor of 200 or more. When the instrument view- 
ed a diffuse homogeneous source of EUV radiation to which it was sensitive, 
less than three per cent of the photons detected reached the photosensitive 
surface by such reflections. 


DATA REDUCTION 


We imposed several constraints on the data selected for analysis and 
discussion here. We required that the line of sight point within 97° of 
the local zenith, which corresponds to the lower edge of the field of view 
being tangent to the Earth's limb. We excluded data taken in or near the 
South Atlantic Anomaly. To avoid contamination by sunlight reflected 

from the spacecraft, we included only data taken when the line of sight 
pointed more than 60° from the Sun. These conditions were met by a total 
of 1105s of data, all taken on July 20 or 21, 1975, from an orbital altitude 


of about 225 km. 
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We binned all data into consecutive 5-s intervals, each, containing 
approximately equal times with the gas absorption cell filled and evacuated. 
The set of 0.1-s telemetered data samples in each bin was filtered to 
eliminate any anomalously high or low count rates, such as might be caused 
by electrical noise or telemetry dropouts. In this process, each 0.1-s 
count-rate sample was examined to see whether it lay within three standard 
deviations of the mean of all other samples in that bin. If not, it was 
discarded. The filtering was repeated until no samples were discarded in 
a complete successive examination of all samples remaining in the bin in 
question. The number of samples eliminated was less than five per cent, 
and was consistent with the number expected from Poisson statistics. 

With the gas absorption cell empty, the count rates detected ranged from 
81 to 3446 counts s 1 . With the cell filled, the corresponding range of 
rates was 12 to 2662 counts s We correlated the reduced cell- full and 
cell— empty data with several parameters upon which airglow intensities 
might depend, including solar zenith angle, line-of-sight elevation, and 
angle between line of sight and the Sun. 

The only correlation we found was with solar zenith angle, which is 
also the angle in arc on the Earth's surface between the subsolar point 
and the spacecraft position. In Figure 2 we show the cell-empty and cell- 
full count rates plotted against solar zenith angle. Statistical count-rate 
errors axe comparable to the size of the dots plotted. The straight lines 
are linear regression fits for the two cases . This dependence of the 
observed fluxes on spacecraft position is qualitatively consistent with 
the model of He I 584 A scattering from an optically thick neutral helium 
geocorona developed by Meier and Weller (1972) . 
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As shown in Figure 2, the cell-full count rate tracked the cell- 
empty count rate remarkably well. Although the fluxes observed varied 
by an overall factor of twenty, the cell-full count rate remained at about 
0.3 of the cell-empty count rate throughout. 


DISCUSSION AND CONCLUSIONS 


The gas absorption cell was designed to absorb all 584 & radiation 
scattered from the Earth's geocorona, even when that radiation was Doppler- 
shifted by the 8 k s ^ speed of an orbiting spacecraft. A nominal 
thermosphere temperature of 1000°K (Jacchia 1971) corresponds to' thermal- 
Doppler broadening of the scattered 584 & line to 3 m2. FWHM, and an addition- 
al 8 km s ■*" speed along the line of sight corresponds to a further Doppler 
shift of about 16 m2. Thus even if the instrument line of sight pointed 
parallel to the orbital motion, the geocoronal helium line would have been 
completely absorbed. For most of the data discussed here, the component of 
orbital velocity along the line of sight was substantially less than 8 km s \ 
The temperatures necessary to thermal-Doppler broaden the geocoronal 584 2 
line to a FWHM comparable to the 99 per cent absorption width of the gas cell, 
are at least an order of magnitude greater than the 1000°K heretofore assumed. 
We are therefore confident that the helium glow detector gas absorption cell 
did in fact absorb all geocoronal 584 2 radiation. 

Radiation at wavelength below the helium continuum edge at 504 2 would 
also be in part absorbed by the gas cell, with the absorption cross-section 
varying in proportion to wavelength cubed. We do not believe such radiation 
could have contributed appreciably to either the cell-full or cell-empty fluxes. 
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We are not aware of any reports of dayglow measurements from 300 to 504 A 

other than broadband photometer studies. However, the solar EUV spectrum 

(Timothy 1976, Hinteregger 1976) shows only one strong line in that range 

likely to be scattered resonantly by atmospheric species. That line is 

© 

the ground state He II transition at 304 A. 

o 

Kumar, Bowyer and Lamp ton (1973) measured 304 A dayglow flux with a 
photometer from altitudes of 120 to 190 km, at a solar zenith angle of 
42°, and after correction for atmospheric absorption found an overhead 
flux of 15 rayleighs. This flux would cause our instrument to count at 
less than 10 counts s \ which is negligible compared to nearly all of our 
data. 

Schmidtke (1976) gives a total solar He I continuum flux comparable 

to the 584 X line emission: that is, about 10 photons cm s . Standard 

14 15 

atmospheric models (Jacchia 1971) indicate approximately 10 - 10 

neutral helium atoms per square centimeter above 220 km. Hence, assuming 

—18 2 

a nominal absorption cross-section of 2 x 10 cm , the instrument would 
observe 0.2-2 rayleigh of atmospheric emission near the neutral helium 
continuum edge, caused by helium photoionization due to the solar continuum 
and subsequent recombination. 

As shown in Figure 1, the instrument sensitivity falls rapidly below 
550 &, so that two rayleighs of flux at 500 X. could cause at most a few 
counts per second. The maximum observed cell-full flux, as shown in Figure 
2, was more than two orders of magnitude greater. 

Both theoretical work (Freeman and Paresce 1978, and references 
discussed therein) and experimental results (Freeman et al. 1978, and 
references therein) indicate that the total solar 584 & flux scattered from 

the interstellar medium passing through the solar system is no more than a 
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few rayleighs for the view directions under consideration here. This 
flux is also much too small to explain our observations. . 

The preceding analysis leads to three conclusions. First, by 
subtracting cell-full from cell-enpty count rates, we conclude that 
the daytime He I EUV geocoronal flux intensity spans the range 2-70 
rayleighs. This value compares favorably with the 25 rayleighs of 584 A 
flux reported by Christensen (1976) between 210 and 240 km at a solar 
zenith angle of 74°. If the 304 A line 'is in fact the -only strong day- 

O 

glow line from 300 to 504 A, then our measurement is not contaminated 
by any radiation whose wavelength differs from He I ground-state lines by 

O 

more than a few mA. 

Second, inasmuch as the gas cell was designed to absorb all geocoronal 
584 £ flux, the substantial 500-700 £ flux observed through the filled cell 
must have some other origin. Interplanetary 584 £ flux might be Doppler- 
shifted sufficiently to be transmitted by the gas cell, but the theoretical 
work cited indicates that this flux is much too weak to explain our cell- 
full data. We cannot calculate an absolute intensity for this radiation, 
since its wavelength or wavelengths, and hence the correct calibration factors 
for our detector, are not known. 

Third, the strong correlation between the cell-full flux and such a 
local parameter as spacecraft position, suggests that this flux originated 
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near the Earth. The intensity of this radiation varies with the observer’s 
position in a manner similar to that of the 584 & flux, at least for the 
range of observation positions and directions considered here. 

Kumar, Bowyer and Lampton (1973) have' reported detection of dayglow 
flux other than at 584 8, within the band 500-800 8. Carruthers and Page 
(1972, 1976) and Christensen (1976) have reported possible detections of 
dayglow spectral lines between 600 and 700 8 with intensities comparable to 
He I 584 8.. Delaboundiniere (1977) has calculated . theoretical emission rates 
for 0 II dayglow lines at 537 and 581 8, at a solar zenith angle of 0°, and 
predicts fluxes of 100-200 rayleighs at the altitudes from which we observed. 
Any of these results might explain our cell-full fluxes. More detailed 
investigations with high-resolution spectrometers will be needed to define 
conclusively the spectral distribution of the dayglow radiation in this band. 
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FIGURE CAPTIONS 

Figure 1: Measured end-to-end photometric efficiency of photometer, 

extrapolated between data points, with typical error bars. 
Figure 2: Cell- full (lower) and cell-empty (upper) count rates versus 

solar zenith angle. Calculated linear regressions fits shown 
by diagonal lines. Typical error bars are given. Five 
points above 2200 counts sec ^ are not shown. 
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ABSTRACT 


We have developed a model for resonant scattering of solar extreme 
ultraviolet light from atoms of the local interstellar helium passing 
through the solar system, which is valid even if the interstellar gas 
is hot. The model' features numerical integration of the phase-space 
number density function of Dariby and Camm (1975) . We compare this model 
with the modified cold-gas model .widely used previously, for values 
of the model parameters presently believed to be nominally correct for 

O 

5S4 A scattering from neutral helium. The new model predicts significantly 
different observed patterns of scattered flux than the old one. The 
new model predicts substantial variations of observed flux with temperature, 
in circumstances in which the old model predicts no variation at all. The 
modified cold-gas model is invalid for analysis of helium observations 
whose field of view lies’ near the downstream symmetry axis of the gas 
motion, and for the determination of the temperature of the distant gas. 

It may be of limited use for other observations, for analyses which depend 
only on the relative intensities of flux observed from different parts of 
the sky. Much previous work must be reexamined' in light of these new 


constraints. 
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INTROPUCTION 

The nearby interstellar medium (ISM) may be directly observed from 
a rocket or spacecraft by solar extreme ultraviolet (EUV) light resonantly- 
scattered from its constituent gas atoms. During the last fifteen years, 

‘a number of authors (see Fahr 1974, Bertaux _et al . 1976, and Freeman el al . 
1976a, and references cited therein) have made such observations, either 
of hydrogen at 1216 A or of neutral helium at 584 A, and have attempted to 
deduce from them properties of the nearby ISM. Several of these quantities, 
notably temperature and number density, are difficult to obtain by other 
means , 

There are no closed formulae whereby these properties may be directly 
calculated by simple substitution of the observed data, All analysis has 
been by' model fitting, In these models, the intensity I of resonantly 
scattered flux which is seen' by an observer with a specified line of sight 
is 'expressed' as an integral: 

(1) I = /n(v,s) F(v a ,s) a(v A , v, p(6(s).) d^vdv^dv.jds 

In this equation, s is path length along the line of sight, v is gas atom 

velocity at the point s,' v is. photon frequency in an inertial frame 

-* 

attached to the Sun, and v is the photon frequency in an inertial frame 

i 

attached to the observer. The function n(v,s) is the phase-space number 

density of scattering atoms at s, a(v , v, v.) is the cross-section for 

« i 

atoms of velocity "v to scatter solar photons of frequency toward the 
observer with frequency 0(s) is the angle between the Sun, the point s, 
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and the observer, p(0(s)) is a phase, .function, and F(v A ,s) is- the solar 

EUV intensity at s. In practice, equation (1) must also include the 

sensitivity of the observing instrument as a factor, and if that 

sensitivity varies significantly over the breadth of the resonant line 

in question, it must appear within the integral, Integrals of this form 

have been described in greater detail in Fahr (1974), Bertaux gt al , (1976) ? 

Freeman j=£ al, (1976a), and in references cited therein, 

Within equation (1) , the forms of a and p are .well-known from atomic 

physics, and F is generally described by a Gaussian of given FWHM and 

intensity, but the proper evaluation of the phase-space number density 

n(v,s) has posed substantial theoretical and numerical difficulties. This 

function has often been replaced by an approximation derived from the 

so-called cold model, in which the thermal speeds of individual gas atoms 

are assumed to be small compared to the bulk' speed V of the distant ISM, 

' B 

When a first-order correction for thermal sppeds which are not small 
compared to the bulk speed is applied to the region downstream from the Sun 
with respect to the approaching ISM (Feldman, Goldstein and Sherb 1972; 

Blum, Pfleiderer and Wulf-Mathies 1975), the model is called the modified 
cold model. It has been used to analyze most data so far published. We 
shall now demonstrate that for analysis of helium observations, the 
modified cold model is wrong by significant amounts in many cases of interest. 


THE MODIFIED COLD MODEL 

The cold model ,phase r space density function derives from the assumption 
.that the distant ISM has a temperature of zero; that is, that all atoms 
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move iii the same direction with the same spaed. Given a- point near the 
Sun, not' located on the axis of symmetry, it can be shown that only two 
orbits with the, zero-temperature initial conditions cross that point. 

The number density function is as given, for example, by Blum, P.fleiderer 
and Wulf-Mathies (1975), equations (4) and (17), 

If the selected- point is directly upstream of the. Sun only one orbit 
crosses it, but in the zero-temperature approximation ' every orbit crosses 
the downstream symmetry axis . This latter condition leads to an unphysical 
density spike at the downstream axis, and the cold model thereby falsely 
predicts a bright, localized EUV flux scattered from that region, Feldman, 
Goldstein and Sherb (1972) and Blum, Pfleiderer and Wulf-Mathies (1975) 
have published formulae for altering the number density function close to 
the downstream axis to take .first-order account of non-zero atomic thermal 
-speeds. With the -addition of such formulae, the cold model becomes the 
modified cold model. 

In this model, the three-dimensional integral over gas velocity in 
equation (1) collapses to a two-term sum, and Dirac delta functions in 
a(v ft , v, v ± ) allow the integrations on v* and v. to be performed immediately, 
The evaluation of the flux I in this model is thereby reduced to a one- 
dimensional numerical integration along the line- of sight, which can be 
accomplished quickly on even a small computer. 

A HOT MODEL 

.The distant ISM is better described, however, by a Maxwell-Boltzmann 
•distribution of temperature T, moving with bulk speed in a certain direction. 
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This initial condition, makes the analytical description of the phase-space 

i «v •* 

number density function n(v,s) a much more difficult problem, Furthermore, 
the use of such numerical techniques as Monte Carlo methods to evaluate 

requires large amounts of time on even the fastest of present-day 
computers. It is therefore not surprising that most workers in the study 
of the local ISM continued to use the modified cold model, even when early 
results indicated that thermal speeds might in fact not be small, 

In 1975 Wallis pointed out that an analytic solution for n(v ? s) under 
the initial conditions just given had been obtained by Danby and Camm 
in 1957, in the context of a different astrophys'ical problem, The Danby 

A 

and Camm expression for n(y,s) contains an exponential of a complicated 
trigonometric function which’ is difficult to manipulate analytically, 
Furthermore, although Dirac delta functions in a(v*, y, V.) allow two of 
the integrations - in equation (1) to be performed immediately, the. calculation 
of an observed flux using the Danby and Camm n(v,s) nevertheless requires 
a four-dimensional integration, three dimensions more than with the modified 
cold model. 


We have developed a program to evaluate equation (1) by numerical 

integration, using Che Danby and Camm number density function, modified 

2 

by the inclusion of the factor expX-g r 0/V p) which was used by preceding 

° E B ’ 

authors (e.g, Blum, Pfleiderer and Wulf-Mathies 1975)- to allow for the loss 


of gas atoms by ionization due to solar EUV,- In this factor, g is the 

? o 

rate at which gas atoms are ionized at one astronomical unit r_ from the 

E ' 

Sun, 0 is the angle through which a gas atom has orbited about the Sun, and 
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p is the impact parameter of the orbit in question. We refer to this model 
as. the hot model. ' Meier has developed a similar model, also using the 
Danby and Camm number density, and has published a description of some 
of the features of its predictions (Meier 1977), Unfortunately, his 
published results so far include only a small number of observer positions 
and lines of sight, and make no direct comparison with the cold model, 

Our hot model vises Simpson's rule to perform all four integrations.. In 
testing the program, we varied the number and spacing of points at which the 
integrand was evaluated, until we observed an apparent asymptotic limit 
to the calculated value of the flux I, For all further work, we used a 
number and spacing of points which produced a value of I within ten per cent 
' of the observed limit . 

We checked the program validity by using several limiting cases' 
singly and in combination. At T = 0, the hot model predictions should equal 
the modified cold model predictions. Since both our computer codes use 
the expression 1/T repeatedly, we could not simply ' substitute T = 0. in the 
program. Instead, we set all other parameters at nominal values and compared 

the predictions of our two models as T- was reduced by Increments to 1 values as 

/ 

low as 300 °K. At these low temperatures, for most lines of sight, the 
models agreed to within about ten per cent, which is the empirically 
determined accuracy of the numerical integration routines used. For lines 
of sight close to the downwind symmetry axis the disagreement at low temperature) 
was much greater. Inasmuch as the analytic forms of both' models predict 
divergent atomic number densities on the downwind axis as T approaches zero, 
this result is .not surprising. Our programs were designed with substantially 
higher temperatures in mind. 
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We also evaluated the hot model in several other limiting cases in 
which some of its features can be checked against analytic expressions. For 
example, with solar mass and radiation pressure set equal to zero, and -with 

O 

solar 584 A line width very large compared to the Doppler shift caused- by 
helium thermal speeds, the hot model predicted an observed line whose 
Gaussian profile matched that calculated - from the Maxwell-Boltzmann distribution 
of the distant gas. The center of this line was . Doppler-shifted by the 
radial velocity between' the gas and the observer. 

Calculations using the hot model also displayed the various symmetries 
imposed by the geometrical symmetries of the physical problem. 

The program provides the total observed flux and the observed line 
profile. It has provision for also calculating the flux and line profile 
observed when a helium gas absorption cell (Freeman et al . 1976b $ Bowyer 
et al. 1977; Fahr et‘ al^ 1977) is interposed between the observer and the 
line of sight. We have so far used the program only to investigate the 
behavior of neutral helium as observed via the resonant scattering of solar 
EUV at 584 A. Each evaluation of the integral (1) requires three to fifteen 
minutes of central processor time on a CDC 6400 computer, 

•MODEL PARAMETERS 

— — r- v •i v r'T V, K ■ 1 » 

Fourteen physical parameters determine the predictions of the 
models. These are three parameters describing the bulk velocity of the 
distant gas, the distant gas temperature T and neutral helium number 
density n , the rate 0 at which the gas is ionized at one A.U, from the 

eo q 

Sun, 'the FWHM AA of the solar resonance line, the solar line flux F at one - 
A.U. from the Sun, the- ratio p of the pressure of light to the force of 
solar gravity oh a gas atom, three parameters ’characterizing the position 
of the observer and two describing the line of sight (LOS) of the observation, 



-7- 


In order to reduce the amount of calculation required to compare the models, 
we restricted several of these parameters to values or‘ narrow ranges near the 
values presently believed corrector to values commonly used in previous 
work. 

We assumed that the distant ISM moves toward the direction a = 72°, 

S = +15°, first given by Weller and Meier (1974) and later confirmed by 
Bertaux et al. (1976). Both models are exactly linear in both and F , 
so that the exact values selected for these parameters alter the predictions 

Q 

only by an overall scale factor. We selected n^ = 0.01 cm -J and F Q = 10 
—7 —1 

cm • s , which are numerically convenient values near the ranges recently 
reported, (Paresce, Bowyer and Kumar 1973, 1974; Weller and Meier 1974; 

Freeman et al, 1977) . We assumed = 6,8 x 10”^ s - ^ and AX = 0,12 A, 

These values are within the ranges reported in recent studies of solar 
properties (Delaboudiniere and Crifo 1975; Doschek, Behring and Feldman 1974; 
Milkey, Heasly and Beebe 1973; Banks and Kockarts 1973), We assumed 
y = 0 throughout. We considered only observers along the Earth's orbit, 
and in general constrained the LOS to the celestial hemisphere opposite 
the Sun. 

RESULTS 

O 

Downwind Sky Map, In Figures 1 and 2 we show maps of 584 A skyglow . 
predicted respectively by the modified cold model and .hot model ? for an 
observer at. the point on the Earth's orbit that is closest (7 degrees of 
arc) to the downwind symmetry axis of the gas motion, for distant gas bulk 
speed 20 km*s -3 - and temperature 10 4 K, Several different kinds of observations 
(Bertaux et al, 1976; Freeman et al. 1976a; Adams and Frisch 1977) indicate 
that is in fact within a few km»s“^ of 20. A temperature of 10^ K is the 
approximate mid-point of the range of gas temperatures determined by recent 
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authors (Paresce, Bowyer and Kumar 1973, 1974; Weller and Meier 1974; 

Bertaux et al«. 1976), The maps are symmetric about the line a = 72° , The 
isophotes for 584 A flux level's in rayleighs are interpolated between 
specific predictions spaced at most 22.5 degrees of arc apart. 

■ Both models predict a pronounced maximum in intensity near the downstream 
direction a = 72 ° $ = +15°, due to the focusing effect of solar gravity. 

The maximum predicted by the modified cold model is more sharply defined than 
the one predicted by the hot model. Furthermore, the modified cold model 
predictions plotted span a greater range than the hot model predictions. 

On the modified cold model map, the highest flux predicted is slightly more 
than twice as great as the lowest. The corresponding factor for the. hot 
model map is 1.5, 

The hot model map shows systematically higher fluxes than the modified 

cold model map, by a factor of from 1,3 to 2.0, Such differences are 

plausible. For example, the' effect of ionization on the predicted flux 

depends on the shapes of the orbits whereby the scattering atoms reached 

the line of sight, and -the effect of the solar line width depends on the 

radial component of gas ‘orbital motion. At T = 10^ K and Yg = 20 km-s"!, 

1 / 2 * 

a helium atom far from the Sun moving with transverse speed V = (2kT/m) 

- 6 km*s _1 has .an orbit vastly different from the one assumed by the 
modified cold model. Such thermally dispersed orbits are 'properly treated 
by the hot model but not by the modified cold model, 

* * O 

Crosswind Sky Map . In Figures 3 and 4 we show maps of 584 A glow 
predicted by the modified cold model and hot model respectively, for an 
observer at a point on the Earth's orbit that is at 90° to the symmetry 
axis. As before, the distant gas bulk speed is 20 km^s" 1 and the temperature 
is 10^ K. The predictions are symmetric about a great circle containing 
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the antisolar point and the direction of gas motion, but the rectangular 
coordinates used in the figures obscure the symmetry. 

Both models predict a localized region of high flux between the down- 
stream- direction of gas motion and the Sun, as is expected from the focusing 
of the gas orbits. However, as in the downstream 'case, the modified cold 
model predicts a markedly sharper region of maximum intensity than does the 
hot model. The predictions of the two models agree within fifteen per cent 
at the point of maximum intensity, but as the line of sight moves away 
from that point they diverge. At the north and south celestial pole^ 
and in' the direction from which the gas approaches, the hot model predicts 
about three times as’ much flux as does the modified cold model, 

The -Effect of Temperature, I-n order to investigate the effect of 

temperature on- the predictions of the two models, we restricted the line of 

sight to directions pointing approximately radially outwapd from the Sun, 

We considered observers located" at five points on. the Earth ^s orbit: as 

/ 

close to the upwind and . downwind symmetry axes. as possible, at the cross wind 
position- discussed above, and at approximately 135° and 160° from the 
upwind symmetry axis. With the LOS pointing approximately outward from the 
Sun from each observer position, we evaluated the model predictions while 
varying T from 625 to 80,000 K, We fixed Vg at 20 kmes"^ and set all other 
parameters to the values discussed above. We show the results 0 £ these 
calculations in Figure 5 for the modified cold model and in Figure 6 for the 
hot model. 

The features of the modified cold model calculations are by now well 
understood. Temperature affects its predictions .only in a region close to 
the downstream symmetry axis, whose exact size, increases with T in the manner 
given by paraboloidally interpolating to the density given by Blum, Pfleiderer 



- 10 - 


and Wulf-Mathies (1975) equation (A17) . The downstream LOS always lies 
within this region, hence the corresponding predictions vary with temperature, 
The high flux at low T indicates sharp focusing. The region where temperature 
variations occur does not expand to include the 160 degree LOS until T is 
between 10^ and 2 x 10^ K: thus predictions at 160 degrees are constant for 

lower temperatures. All other lines of sight never penetrate the temperature- 
sensitive region, hence the model predictions remain constant as T changes, 

In contrast, the hot model exhibits considerable variation of 
predictions with temperature for all lines of sight. Note that even as 
far as 90 degrees from the downwind direction, varying T over the 1500 
to 20,000 K range produces a fifty per cent change in predicted flux. 

In the -downwind direction at low temperature, the predictions of the 
two models differ -substantially , We advise caution In ascribing great 
significance to this difference, because both models encounter unphysical 
infinities in the limit T=0. In the modified cold models, pa rticl e number 
density on the downwind axis diverges, and in the hot model, phase-space 

number density approaches a sum of Dirac delta functions. Inasmuch as our 

* * » . * * 

coding . for both of these models was designed with temperatures above 

3000 K in mind, it is possible that at substantially lower temperatures, 

the coding for one or both models produces inaccurate results, 

IMPLICATIONS FOR RECENT DATA ANALYSES 

After performing the tests and general calculations just discussed ^ 

0 

we used the .hot model to check some recently published 584 A 

observations which have previously been analyzed with the modified cold 
model. We first selected the satellite observations reported by Weller and 
Meier (1974) . These authors presented data gathered during almost a year 

Q 

by a '584 A photometer whose LOS- always maintained a constant 25-degree 
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angle with respect to the plane containing the Sun, the Earth, and the 
north celestial pole. Their instrument revolved rapidly about the- normal 
to that plane, so that its LOS traced a cone of 130° apex angle. The 
cone in turn revolved at one revolution per year. 

The greatest variations in the Weller and Meier data take place 
approximately at declination +10°, as show in Weller and Meier (1974), 

-’Figure 2. We selected values of Weller and Meier coordinate "d", which 
describes the location and orientation of the cone, for which the 584 A 
isophotes in that Figure cross 6 = +10°, and evaluated our hot model for 
the resulting observer positions and lines of sight. We used temperatures 
of 5000, 10,000 and 20,000 K, and set all other parameters as described 
above. 

The results of our calculations are shown in Figure 7, in which the curve o 
short dashes is drawn through the data points taken from Weller and Meier 
and the curves labeled by temperature depict our hot-model predictions, 

In order to display the difference in shape of the curves clearly, we have 
multiplied each of our curves by a scale factor, to produce the same 
maximum intensity as the Weller and Meier data. This multiplication simply 
corresponds to an adjustment of the overall linear factors n and F in the 
hot model. 

Although each calculated curve might be said to have the same general 
shape as the Weller-Meier data, curve, there are two systematic differences, 
between theory and observations. First, the shapes of the curves are 
' distinctly different near the region of maximum flux. The data resemble 
a Gaussian, but the calculated curves are more nearly triangular, with 
pointed tips and substantially broader bases. Second, none of the calculated 
curves display as high a ratio of maximum to minimum flux as do the data, 

We therefore are reluctant to state that the hot model fits that we have 
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obtained favor a specific temperature. Weller and Meier used the modified 
cold model to obtain a range of good fits between 5000 and 20,000 KL, with the 
same values for the speed and direction of gas motion that we have used. 

We conjectured that our calculations might have been done with 
too high a value of the ionization rate $ Q . If the selected value of 
6,8 x 10”® s”^ were greater than that actually present, then the model 
would predict excessively low helium number densities near the downwind 
symmetry axis, and consequently low flux levels, However, a few calculations 
(not shown) made with 6 0 = 0 raised the ratio of maximum to minimum flux 
only half way to that observed by Weller and Meier, 

Schmidtke (1976) has recently reported variations of the solar 584 A 
flux F Q of as much as 30 per cent on a time scale of days to months. If 
had been 30 -per cent higher for the a = 80° data than for the a = 260° 
data, then the theory and observations would be in better agreement. 
Unfortunately , Schmidtke’s data do not include the time period in late 1972 
when the Weller and Meier instrument was observing the region of high 
584 A flux, so that no simple adjustment for varying F q is possible. 

We have not determined what variations of other model parameters, such 
as the bulk speed and direction of motion of the distant gas , will produce 
a best fit of the hot model to the Weller and Meier observations. To do 
so would require a much more extensive analysis, using all of the Weller- 
Meier data instead of the small subset that we have chosen here. But our 
calculations clearly demonstrate that the hot and modified cold models can 

O 

give substantially different results when used to analyze helium 584 A data. 
The differences between the models are therefore of considerable practical 
importance, for the derivation of ISM parameters from observations. 
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Freeman et al . (197 6a, b) used a 584 A photometer equipped with a 

helium gas resonant absorption cell to study radiation scattered from 

a ~ 50°, 6 = 39°. In their analysis, they used the ratio e Q of cell-full 

to cell-empty flux as a parameter for model fitting, and thereby reached 

conclusions independent of the exact values of n m and F . They obtained 

a lower limit, at a 95 per cent confidence level, of 10-15 km* s"^ on the 

distant gas- bulk speed V . We have used the hot model in an attempt to fit 

b 

their data, with T = 5000, 10,000 and 20,000 K and Vg = 5, 10, 15, 20, 30, 
and 40 km»s“^. All other parameters were set at the yalues discussed 
above. In Figure 8 we show how the hot model calculations fit the data 
of Freeman si al. . The horizontal dashed line is the observed value of 
£ q , and the shaded area depicts two standard deviation limits thereon. 

The curves labeled by temperature are the model predictions, which 

vary by up to 30 .per cent over the range 5000-20,000 K, The lower limit 
on Vg indicated by the hot model, at a- 95 per cent confidence level, is 
17-20 km^s* - "**, which is in better agreement with the recent results of 
Adams and Frisch (1977) and Bertaux et. alj (1976) than was the limit derived 
using the modified cold model, 

CONCLUSIONS 

Figures 1-4 demonstrate that when the modified cold model and the 
hot model predictions for helium atoms are evaluated for physical 
parameters in the range presently believed correct, the hot model predicts 
broader intensity maxima and smaller total variations In flux level than 
does the modified cold model, The shapes of the intensity maxima are 
sufficiently complicated that the difference in shape is difficult to 
quantify* but the calculations presented here indicate that the patterns 

O 

of 584 A glow predicted by the two models are markedly different within 
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60 to 90° away from the projection of the downwind symmetry axis on the 
celestial sphere. We conclude that the modified cold model should not be 
used to analyze observations of this region. 

Figures 3 and 4 indicate that at greater distances from the downwind axis, 

the shapes and spacings of the isophotes predicted by the two models are 

more similar, although the hot model predicts several times as much flux 

as the modified cold model. The modified cold model may possibly be of 

use for interpreting' observations of this part of the sky, provided the 

analysis used depends only on the relative intensities of flux observed from 

different points. Future observers of 584 & helium glow should address 

the question of model validity more specifically, 

* 

Figures 5 and 6 demonstrate that the hot model predicts substantial 
variations of observed flux with temperature, on lines of sight for which 
the modified cold model predicts no variation at all. Furthermore, even 
when the modified cold model does predict a variation of flux with 
temperature, the hot model predicts variations of a different kind. 

The use of the modified cold model to determine the precise temperature 
of the local ISM therefore appears invalid, * 

O 

The conclusions drawn by previous observers of 584 A flux who have 
used the modified cold model to interpret their data must be reassessed , 

The hot model does not reproduce the fit that Weller and TIeier (19 74) 
obtained using the modified cold model. We suggest that further 
study is in order, to determine the cause of the discrepancy 
and the actual values of the' ISM parameters that the data imply-, 

Foint-by-point calculations of the particle number densities and velocity 
dispersions predicted by the two models would be useful in this context. 

When the hot model is used to analyze the data of Freeman et al. (1976) the 
derived lower limit to gas bulk speed is sJLightly nearer to the values of 
obtained by other means than was the case with the modified cold model. 
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The 'theoretical question of the behavior of the local ISM is far from 
closed. Both the modified cold model and the hot model derive from two 
assumptions. First, the flow of the ISM through the solar system is 
assumed to be collisionless, so that each atom proceeds along a hyperbolic 
orbit determined by celestial mechanics. Second, the- nearby ISM is assumed 
sufficiently tenous that only singly-scattered' photons will reach the 
observer, Wallis (1973, 1974, 1975) has questioned both of these 
assumptions. Further analytic and numerical study of their validity may 
be necessary before experimenters can have complete confidence in the 

' O 

interpretation of 584 A data. 
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FIGURE CAPTIONS 
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Figure 4 : 


Downwind Map-Modified Cold Model 

- O 

Modified cold model predictions of interplanetary 584 A glow 
seen from the downwind side of the Earth's orbit, for V = 20 km s 

E> 

4 

T - 10 K. Isophotes in rayleighs interpolated between calculated 
predictions. The map is symmetric about a = 72°. Direction 
toward which gas flows indicated by 0. 

Downwind Map-Hot Model 

O 

Hot model predictions of interplanetary 584 A glow seen from 

the downwind side of the Earth’s orbit, for = 20 km s \ 

i5 

4 

T ~ 10 K. Isophotes in rayleighs interpolated between calculated 
predictions. The map is symmetric about a = 72°. Direction 
toward which gas flows indicated by 0. 

Crosswind Map-Modified Cold Model 

O 

Modified cold model predictions of interplanetary 584 A glow 
seen from the Earth's orbit midway between upwind and downwind 
symmetry axes, for V„ = 20 km s 1 , T = 10^ K. Isophotes in 

-D 

rayleighs interpolated between calculated predictions . Gas 
flows from direction indicated by ©, toward direction indicated 
by 0. Antisolar point indicated by &. 

Crosswind Map-Hot Model 

0 

Hot model predictions of interplanetary 584 A glow seen from 

the Earth's orbit midway between upwind and downwind symmetry. 

-1 4 

axes, for Vg = 20 km s , T = 10 K. Isophotes in rayleighs 
interpolated between calculated predictions . Gas flows from 
direction indicated by €, toward direction indicated by 0. 


Antisolar point indicated by a. 



Figure 5: 


Figure 6 : 


Figure 7: 


Figure 8: 


Effect of Temperature-Modified Cold Model 

O 

Modified cold model predictions of interplanetary 584 A glow 
seen from various points on the Earth's orbit, looking 
radially outward from the Sun, as a function of temperature , 
for V- = 20 km s \ Curves interpolated between calculated 
predictions. Curve labels are angles between line of sight 
and upwind symmetry axis. 

Effect of Temperature-Hot Model 

O 

Hot model predictions of interplanetary 584 A glow seen from 
various points on the Earths's orbit, looking radially outward 
from the Sun, as a function of temperature , for = 20 km s 
Curves interpolated between calculated predictions. Curve 
labels are angles between line of sigfit and upwind symmetry axis . 
Weller-Meier Data Fit 

Use of hot- model to fit data of Weller and Meier (19 74). Line 
of short dashes: interpolation between Weller and Meier data 

points. Other lines: hot model predictions for the stated 

temperatures, scaled to the same maximum flux observed by 
Weller and Meier. The coordinate a describes the Weller and 
Meier viewing geometry. 

Freeman et al. Data .Fit 

Use .of hot model to fit data of Freeman et al. (1976a). 

Dashed line: observed ratio of cell-full to cell-empty flux. 

Shaded area: two standard deviation limits on ratio. Labeled 

curves: hot model predictions for the stated temperatures. 


as a function of distant gas bulk speed. 
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